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Abstract. Astrophysical fluids are generally turbulent and 
this preexisting turbulence must be taken into account for 
the models of magnetic reconnection which are attepmted 
to be applied to astrophysical, solar or heliospheric environ- 
ments. In addition, reconnection itself induces turbulence 
which provides an important feedback on the reconnection 
process. In this paper we discuss both theoretical model and 
numerical evidence that magnetic reconnection gets fast in 
the approximation of resistive MHD. We consider the rela- 
tion between the Lazarian & Vishniac turbulent reconnec- 
tion theory and Lapenta's numerical experiments testifying 
of the spontaneous onset of turbulent reconnection in sys- 
tems which are initially laminar. 
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1 Introduction 

Astrophysical plasmas are known to be magnetized and tur- 
bulent. Magnetization of these fluids most frequently arises 
from the dynamo action to which turbulence is an essen- 
tial component (see Schekochihin et al. 2007). The drivers 
of turbulence, e.g. supernovae explosions in the interstellar 
medium, inject energy at large scales and then the energy 
cascades down to small scales through a hierarchy of ed- 
dies spanning up over the entire inertial range. The famous 
Kolmogorov picture (Kolmogorov 1941) corresponds to hy- 
drodynamic turbulence, but, as we discuss further, a quali- 
tatively similar turbulence also develops in magnetized flu- 
ids/plasmas. The definitive confirmation of turbulence pres- 
ence comes from observations, e.g. observations of electron 
density fluctuations in the interstellar medium, which pro- 
duce a so-called Big Power Law in the Sky (Armstrong et al. 
1994, Chepurnov & Lazarian 2010), with the spectral index 
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coinciding with the Kolmogorov one. A more direct piece 
of evidence comes from the observations of spectral lines. 
Apart from showing non-thermal Doppler broadening, they 
also reveal spectra of supersonic turbulent velocity fluctua- 
tions when analyzed with techniques like Velocity Channel 
Analysis (VCA) of Velocity Coordinate Spectrum (VCS) de- 
veloped (see Lazarian & Pogosyan 2000, 2004, 2006, 2008) 
and applied to the observational data (see Padoan et al. 2004, 
2009, Chepurnov et al. 2010) rather recently. 

Reconnection is a process of changing magnetic flux 
topology and this process has been a challenge to understand 
in highly conductive astrophysical fluids. Naturally, it is im- 
portant to understand the process of magnetic reconnection 
in reallistically turbulent state of the fluid. Observationally, 
it is also known that a high level of turbulence is present as 
outcome of solar flares, which are generally believed to be 
driven by magnetic reconnection. Therefore it is important to 
understand the mutual feedback of turbulence and reconnec- 
tion and provide the connection between the existing theory, 
observations and numerical experiments. This is the goal of 
the present paper. 

We consider both the model of turbulent reconnection sug- 
gested in Lazarian & Vishniac (1999, henceforth LV99) and 
numerical evidence of fast reconnection in a few numerical 
papers by Lapenta and coauthors. We search for the relation 
of the theory and observations. 

The very idea that turbulence can influence magnetic re- 
connection is not new. However, LV99 model is radically 
different from its predecessors which also appealed to the ef- 
fects of turbulence. For instance, unlike Speiser (1970) and 
Jacobson (1984) the model does not appeal to changes of the 
microscopic properties of the plasma. The nearest progenitor 
to LV99 was the work of Matthaeus & Lamkin (1985, 1986, 
henceforth ML), who suggested that magnetic reconnection 
may be fast due to a number of turbulence effects, e.g. multi- 
ple X points and turbulent EMF. However, ML did not realize 
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the key role of played by magnetic field-line wandering] and 
did not obtain a quantitative prediction for the reconnection 
rate, as did LV99. From the numerical simulations in ML 
obtaining the actual reconnection rate is not straightforward. 
They do not discuss the onset of reconnection from the lam- 
inar case, either. On the contrary, we shall demonstrate the 
results of simulations where both of these effects are studied 
and quantified. 

In what follows, we discuss magnetic reconnection in tur- 
bulent fluid in §2, provide the numerical confirmations of the 
predictions of the LV99 model in §3, consider the sponta- 
neous onset of reconnection in MHD simulations in §4 and 
the role of flow pattern in §5. We discuss the implications 
of the expected flares of reconnection in §6 and provide our 
summary in §7. 

2 Reconnection can be Fast in Turbulent Plasmas 

Astrophysical plasmas are often highly ionized and highly 
magnetized (Parker 1970). The evolution of the magnetic 
field in a highly conducting fluid can be described by a sim- 
ple version of the induction equation 

OB 

— =Vx(DxB-t)VxB), (1) 

where B is the magnetic field, v is the velocity field, and rj is 
the resistivity coefficient. Under most circumstances this is 
adequate for discussing the evolution of magnetic field in an 
astrophysical plasma. When the dissipative term on the right 
hand side is small, as is implied by simple dimensional es- 
timates, the magnetic flux through any fluid element is con- 
stant in time and the field topology is an invariant of mo- 
tion. On the other hand, reconnection is observed in the solar 
corona and chromosphere (Innes et al. 1997, Yokoyama & 
Shibata 1995, Masuda et al. 1994, Ciaravella & Raymond 
2008), its presence is required to explain dynamo action in 
stars and galactic disks (Parker 1970, 1993), and the violent 
relaxation of magnetic fields following a change in topol- 
ogy is a prime candidate for the acceleration of high energy 
particles (de Gouveia Dal Pino & Lazarian 2003, henceforth 
GL03, 2005, Lazarian 2005, Drake et al. 2006, Lazarian & 
Opher 2009, Drake et al. 2010) in the universe. Quantitative 
general estimates for the speed of reconnection start with two 
adjacent volumes with different large scale magnetic field di- 
rections (Sweet 1958, Parker 1957). 

LV99 we introduced a model that included the effects of 
magnetic field line wandering (see Figure 1). The model re- 
lies on the nature of three-dimensional magnetic field wan- 
dering in turbulence. This nature is different in three and two 
dimensions, which provides the major difference between the 
LV99 model and the earlier attempts to solve the problem of 
magnetic reconnection appealing to turbulence (Matthaeus 

'instead a discussion other effects, like heating, was presented. 
Those do not play a role in LV99. 
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Fig. 1. Upper plot: Sweet-Parker model of reconnection. The out- 
flow is limited by a thin slot A, which is determined by Ohmic dif- 
fusivity. The other scale is an astrophysical scale t> A. Middle 
plot: Reconnection of weakly stochastic magnetic field according 
to LV99. The model that accounts for the stochasticity of magnetic 
field lines. The outflow is limited by the diffusion of magnetic field 
lines, which depends on field line stochasticity. Low plot: An in- 
dividual small scale reconnection region. The reconnection over 
small patches of magnetic field determines the local reconnection 
rate. The global reconnection rate is substantially larger as many 
independent patches come together. From Lazarian et al. 2004. 



& Lamkin 1985). The effects of compressibility and heat- 
ing which were thought to be important in the earlier stud- 
ies (Matthaeus & Lamkin 1985, 1986) do not play the role 
for the LV99 model either. The model is applicable to any 
weakly perturbed magnetized fluid, irrespectively, of the de- 
gree of plasma being collisional or collisionless (cf. Shay et 
al. 1998). 

Two effects are the most important for understanding of 
the nature of reconnection in LV99. First of all, in three di- 
mensions bundles of magnetic field lines can enter the re- 
connection region and reconnection there independently (see 
Figure 1), which is in contrast to two dimensional picture 
where in Sweet-Parker reconnection the process is artificially 
constrained. Then, the nature of magnetic field stochasticity 
and therefore magnetic field wandering (which determines 
the outflow thickness, as illustrated in Figure 1) is very dif- 
ferent in 2D and the real 3D world (LV99). In other words, 
removing artificial constraints on the dimensionality of the 
reconnection region and the magnetic field being absolutely 
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straight, LV99 explores the real-world astrophysical recon- 
nection. 

Our calculations in LV99 showed that the resulting recon- 
nection rate is limited only by the width of the outflow re- 
gion. This proposal, called "stochastic reconnection", leads 
to reconnection speeds close to the turbulent velocity in the 
fluid. More precisely, assuming isotropically driven turbu- 
lence characterized by an injection scale, I, smaller than the 
current sheet length, we find 

Vreo « ^ (l/L) 1/2 « U turb (l/Lf 2 , (2) 
y A 

where ui is the velocity at the driving scale and u tU rb is the 
velocity of the largest eddies of the strong turbulent cascade. 
Note, that here "strong" means only that the eddies decay 
through nonlinear interactions in an eddy turn over time (see 
more discussion of the LV99). All the motions are weak in 
the sense that the magnetic field lines are only weakly per- 
turbed. 

It is useful to rewrite this in terms of the power injection 
rate P. As the perturbations on the injection scale of turbu- 
lence are assumed to have velocities ui < Va, the turbulence 
is weak at large scales. Therefore, the relation between the 
power and the injection velocities are different from the usual 
Kolmogorov estimate, namely, in the case of the weak turbu- 
lence P~uf/(IV A ) (LV99). Thus we get, 

where I is the length of the turbulent eddies parallel to the 
large scale magnetic field lines as well as the injection scale. 

The reconnection velocity given by equation ^ is ob- 
tained in MHD limit and therefore it does not depend on re- 
sistivity or plasma effects. Therefore LV99 model predicts 
that for sufficiently high level of turbulence collisionless and 
collisional fluids should reconnect at the same rate. 

3 Fast reconnection in Turbulent MHD Simulations 

Here we describe the results of a series of three dimensional 
numerical simulations aimed at adding turbulence to the sim- 
plest reconnection scenario and testing equation We take 
two regions with strongly differing magnetic fields lying next 
to one another. The simulations are periodic in the direction 
of the shared field (the z axis) and are open in the reversed di- 
rection (the x axis). The external gas pressure is uniform and 
the magnetic fields at the top and bottom of the box are taken 
to be the specified external fields plus small perturbations to 
allow for outgoing waves. The grid size in the simulations 
varied from 256x512x256 to 512x1028x512 so that the top 
and bottom of the box are far away from the current sheet 
and the region of driven turbulence around it. At the sides 
of the box where outflow is expected the derivatives of the 
dynamical variables are set to zero. A complete description 



of the numerical methodology can be found in Kowal et al. 
(2009). All our simulations are allowed to evolve for seven 
Alfven crossing times without turbulent forcing. During this 
time they develop the expected Sweet-Parker current sheet 
configuration with slow reconnection. Subsequently we turn 
on isotropic turbulent forcing inside a volume centered in the 
midplane (in the xz plane) of the simulation box and extend- 
ing outwards by a quarter of the box size. The turbulence 
reaches its full amplitude around eight crossing times and is 
stationary thereafter. 

In Figure 2 we see the current density on an xy slice of 
the computational box once the turbulence is well developed. 
As expected, we see that the narrow stationary current sheet 
characteristic of Sweet-Parker reconnection is replaced by a 
chaotic structure, with numerous narrow peaks in the current 
density. Clearly the presence of turbulence has a dramatic 
impact on the structure of the reconnection zone. In addition, 
we see numerous faint features indicating weak reconnection 
between adjacent turbulent eddies. 

The speed of reconnection in three dimensions can be hard 
to define without explicit evaluation of the magnetic field 
topology. However, in this simple case we can define it as 
the rate at which the x component of the magnetic field dis- 
appears. More precisely, we consider a yz slice of the sim- 
ulation, passing through the center of the box. The rate of 
change of the area integral of — B x — is its flux across the 
boundaries of the box minus the rate at which flux is annihi- 
lated through reconnection (see more discussion in Kowal et 
al. 2009) 

ft (J\B x \ dzd y) -f^B^V-W^^LM 

where electric field is E = v x B— r)j, B xext is the abso- 
lute value of B x far from the current sheet and L z is the 
width of the box in the z direction. This follows from the 
induction equation under the assumption that the turbulence 
is weak to lead to local field reversals and that the stresses 
at the boundaries are weak to produce significant field bend- 
ing there. In other words, fields in the x direction are ad- 
vected through the top and bottom of the box, and disappear 
only through reconnection. Since we have assumed periodic 
boundary conditions in the z direction the boundary integral 
on the right hand side is only taken over the top and bottom 
of the box. By design this definition includes contributions to 
the reconnection speed from contracting loops, where Ohmic 
reconnection has occurred elsewhere in the box and \B X \ de- 
creases as the end of a reconnected loop is pulled through 
the plane of integration. It is worth noting that this estimate 
is roughly consistent with simply measuring the average in- 
flux of magnetic field lines through the top and bottom of 
the computational box and equating the mean inflow veloc- 
ity with the reconnection speed. Following equation Q we 
can evaluate the reconnection speed for varying strengths and 
scales of turbulence and varying resistivity. 
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Fig. 2. Left panel: Current intensity and magnetic field configuration during stochastic reconnection. We show a slice through the middle of 
the computational box in the xy plane after twelve dynamical times for a typical run. The shared component of the field is perpendicular to 
the page. The intensity and direction of the magnetic field is represented by the length and direction of the arrows. The color bar gives the 
intensity of the current. The reversal in B x is confined to the vicinity of y=0 but the current sheet is strongly disordered with features that 
extend far from the zone of reversal. Right panel: Representation of the magnetic field in the reconnection zone with textures. 
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Fig. 3. Reconnection speed versus input power for the driven tur- 
bulence. We show the reconnection speed, defined by equation (4) 
plotted against the input power for an injection wavenumber equal 
to 8 (i.e. a wavelength equal to one eighth of the box size) and a 
resistivity v u . The dashed line is a fit to the predicted dependence 
of P 1 / 2 (see eq. (3)). The horizontal line shows the laminar recon- 
nection rates for each of the simulations before the turbulent forcing 
started. Here the uncertainty in the time averages are indicated by 
the size of the symbols and the variances are shown by the error 
bars. 



In Figure|5]we see the results for varying amounts of input 
power, for fixed resistivity and injection scale as well as for 
the case of no turbulence at all. The line drawn through the 
simulation points is for the predicted scaling with the square 
root of the input power. The agreement between equation (|3]l 
and Figure [3] is encouraging but does not address the most 
important aspect of stochastic reconnection, i.e. its insensi- 
tivity to r\. 

In Figure |4] we plot the results for fixed input power and 
scale, while varying the background resistivity. In this case rj 
is taken to be uniform, except near the edges of the compu- 
tational grid where it falls to zero over five grid points. This 
was done to eliminate edge effects for large values of the re- 
sistivity. We see from the Figure |4]that the points for laminar 
reconnection scale as y/fj, the expected scaling for Sweet- 
Parker reconnection. In contrast, the points for reconnection 
in a turbulent medium do not depend on the resistivity at all. 
In summary, we have tested the model of stochastic recon- 
nection in a simple geometry meant to approximate the cir- 
cumstances of generic magnetic reconnection in the universe. 
Our results are consistent with the mechanism described by 
LV99. The implication is that turbulent fluids in the universe 
including the interstellar medium, the convection zones of 
stars, and accretion disks, have reconnection speeds close to 
the local turbulent velocity, regardless of the local value of re- 
sistivity. Magnetic fields in turbulent fluids can change their 
topology on a dynamical time scale. 

In Kowal et al. (2009) we also studied the dependence 
of the reconnection on the anomalous resistivity, which in- 
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Fig. 4. Reconnection speed versus resistivity. We show the recon- 
nection speed plotted against the uniform resistivity of the simula- 
tion for an injection wavenumber of 8 and an injected power of one. 
We include both the laminar reconnection speeds, using the hollow 
symbols, fit to the expected dependence of rj u , and the stochastic re- 
connection speeds, using the filled symbols. As before the symbol 
sizes indicate the uncertainty in the average reconnection speeds 
and the error bars indicate the variance. We included simulations 
with large, B z — 1, and small, B z — 0.1, guide fields. 



creases effective resistivity for high current densities. The 
anomalous resistivity can be used as a proxy for plasma ef- 
fects, e.g. collisionless effects in reconnection. While it en- 
hances the local speed of individual reconnection events re- 
sults in Kowal et al. (2009) testify that the total reconnection 
rate does not change. In addition, the study of numerical 
effects is presented in Lazarian et al. (2011). Numerical re- 
sistivity decreases with increase of the numerical resolution. 
Therefore, if the numerical effects were influencing the re- 
connection rate, the increase of the resolution would decrease 
the reconnection rates. This is not what is seen in simulations 
(Lazarian et al. 2011). If anything, the reconnection rates 
slightly increase with the increase of resolution which is due 
to the fact that the turbulence proceeds to smaller scales. As 
a result, field wandering over smaller scales is available and 
therefore the outflow region gets slightly thicker as the nu- 
merical viscosity decreases. 

Finally, it is important to give a few words in relation to 
our turbulence driving. We drive our turbulence solenoidally 
to minimize the effects of compression, which does not play 
a role in LV99 model. The turbulence driven in the volume 
around the reconnection layer corresponds to the case of as- 
trophysical turbulence, which is also volume -driven. On the 
contrary, the case of the turbulence driven at the box bound- 
aries would produce spatially inhomogeneous imbalanced 
turbulence for which we do not have analytical predictions 
(see discussion of such turbulence in Beresnyak & Lazarian 
2009). We stress, that it is not the shear size of our numerical 
simulations, but the correspondence of the observed scalings 
to those predicted in LV99 that allows us to claim that we 



proved that the 3D reconnection is fast in the presence of tur- 
bulence. 



4 Spontaneous onset of Turbulent Reconnection in 
Laminar System 

Macroscopically laminar systems can spontaneously transi- 
tion toward unsteady regimes characterised by faster rates 
and presenting features that are closely related to the turbu- 
lent reconnection regimes described above. 

The long recognized natural pathway for reconnection in 



laminar systems is the Sweet-Parker (SP) regime (Sweet 
|1958| |Parker| |1963| >. In SP reconnection a thin and elon- 
gated current layer forms separating two areas of opposing 
magnetic polarity. The plasma and the field flows towards 
the current layer in the direction normal to the elongated side 
and outflows along the axis of the current. Figure Qja shows 
the typical flow pattern. 

The SP current layer, like any other current layer, is in it- 
self unstable to the tearing mode (Biskamp 1993 2000} pro- 
ducing secondary islands. The presence of the flow pattern 
and the finite length of the layer prevent a direct application 
of the textbook analysis of the tearing mode. The first studies 
of the instability of SP layers were reported by Bulanov et al. 



dT979j) and are summarised in recent textbooks (Biskamp 
|1993||200"0] ). The conclusion was that for aspect ratios of the 
SP layer (length of the current channel, A over its width 5) 
exceeding approximately 100, the SP layer becomes unsta- 
ble. The aspect ratio of a SP layer is directly determined by 
the resistivity in the system. Using the parametrization of 
the resistivity with the Lundquist number S, the aspect ratio 



is simply A/S = \/S |Priest and Forbes 



2000). It follows 



that given a sufficiently low resistivity, corresponding to a 
Lundquist number in excess of S = 10 4 , a SP will become 
naturally unstable to formation of secondary islands. 

Low resistivity was not easily accessible to simulations un- 
til recently. At low resolutions, the numerical resistivity ex- 
ceeds the intended low resistivity preventing the exploration 
of truly low resistive regimes. However, in recent years truly 
low resistivity regimes became accessible thanks to modern 
numerical methods and state of the art computational facil- 
ities. Great progress has been made in understanding this 
spontaneous transition from laminar SP reconnection to an 
unsteady regime characterised by a random onset of sec- 
ondary islands. 

The configuration of a SP layer is intrinsically 2D and re- 
quires a 2D analysis ( |Ni et al. |2010[ ). Furthermore the SP 
layer is embedded in macroscopic configurations that further 
affect the stability and evolution (Schindler and Birn 1999). 
However, a reduced ID analysis ( |Loureiro et al.| |2007| ) can 
still take into account the effects of the flow present in the 
SP layer. The flow has a stabilising effect ( |Bulanov et al.| 



1978 ) that prevents the onset for aspect ratios below approx- 



imately 100, but at higher aspect ratios (corresponding to 
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higher Lundquist numbers) the instability is present. The 
growth rate of the secondary instability is found to scale 
as 7 « S 1 / 4 and the fastest growing mode corresponds to a 
number of islands that scales as S 3 / 8 ( Huang and Bhattachar- 
[jee 2010 1. The power is positive: the growth rate of the sec- 



entry of the reconnection region VA,in- F° r a reconnecting 



ondary instability is faster for lower resistivity. Additionally 
as the resistivity is lowered and the aspect ratio is increased, 
the number of magnetic islands increases. The emerging pic- 
ture from linear theory is then clear: as the Lundquist number 
is increased the SP layer becomes progressively more elon- 
gated, with a higher aspect ratio. The more elongated the SP 
layer becomes, the larger the number of secondary islands 
becomes, and the faster they grow after forming. A clear 
indication of transition towards a turbulent regime like that 
covered in the previous sections. 

The prediction of the linear theory with respect to the num- 
ber and growth rate of secondary islands have been tested and 



verified in carefully designed simulations (Samtaney et al 



20091 |Cassak etaL] [2009 Bha ttacharjee et al.||2009| ) . 

Similarly, the threshold for onset (S approximately 10 4 ) 
has been also confirmed (Skender and Lapenta 2010). Addi- 
tionally, the simulation studies have proven that turbulence, 
if present, has a direct impact on the process. The pres- 
ence of pre-existing fluid turbulence changes the onset of 
the secondary island instability (S kender and Lapenta||2010"[ 
|Loureiro et al.| |2009[ ). A SP layer still forms also in pres- 
ence of moderate levels of turbulence but it is more prone 
to become unstable. The stronger seed reduces the thresh- 
old for the onset and the system transitions more quickly to 
an unsteady regime of reconnection characterized by multi- 
ple islands and reconnection sites. Indeed in all simulations 
the instability of the SP layer is seeded by numerical noise. 
Spectral methods avoid such noise completely and inedeed 
the instability cannot develop ( |Ng and Ragunathan 201 1) , 
unless a seed in the form of a low level of turbulence is added 
explicitly. 

An important consequence of the formation of secondary 
islands and the onset of unsteady reconnection is the accel- 
eration of the reconnection process (Lapenta, 2008). Recon- 
nection is usually defined as fast when two conditions are 
met. 

First, the reconnection rate should be independent of the 
mechanism that allows the decoupling of the plasma (and 
specifically the electrons in the case of advanced two fluid 
and kinetic models). In resistive MHD, the mechanisms al- 
lowing such decoupling are resistivity and, possibly, viscos- 
ity. With the onset of secondary islands, simulations observe 
not only an increased rate of reconnection but careful stud- 
ies varying the resistivity in the system have shown that the 
rate becomes insensitive or even independent of resistivity 
altogether (|Lapenta 2008[ |Loureiro et aL 2009 Huang and 
|Bhattacharjee||2010||Uzdensky et al.||2010) . 

Second, for reconnection to be fast in absolute terms, its 
rate must be a significant fraction of the local Alfven speed 
measured with the density and magnetic field strength at the 



Harris sheet, [Skender and Lape nta (2010) report a rate of ap- 
proximately A^f /vA,i n Bi n ~ .015 with peaks of twice that 
value. Note that the value reported by Skender and Lapenta 



(20101 is computed with a different normalisation Alfven 
speed based on the peak rather than the inflow density that 
is a factor of 10 smaller. Here we prefer to translate the re- 
sult in terms of to the actual Alfven speed of the inflowing 
plasma for better comparison with the other litarature. For 
force free equlibria, the peak reconnection rate also nearly 



reaches /vA in B in w .015 (Skender and Lapenta 



2010). 



Similarly, in the case of the reconnection induced by the co- 
alescence of two magnetic islands, Huang and Bhattachar- 



jee 



(2010 1 report a reconnection rate of A^/vaBq sa .01, in 



agreement with the case of the Harris sheet. The reconnec- 
tion rate of unsteady MHD resistivity is therefore somewhat 
lower than that typical of the fastest kinetic reg imes (peak 
rate of AW /vA,i n Bi n ~ -1 (Bim and et al. 2001 1) but is still 
very strong. 

Many of the properties of the fast reconnection process 
following the onset of the secondary islands instability can 
be understood with a simplified approach based on the prop- 
erties of the tearing mode (Bhattacharje e et aL]|2009 1. For an 
initial current sheet with a magnetic field profile 



B(z) = _B tanh(z/a)x 



(5) 



the linear theory based on reduced MHD with uniform re- 



sistivity r) predicts a maximum growth rate ( Schindler 2006 
p.242): 



Imaxa/vA = 0.623 -(n av A /v) 



-1/2 



corresonding to the wave number: 



kmaxd = 1.358- {noavA/v) 1/4 



(6) 



(7) 



In the case of the stability of a SP layer, one can simply 
use the same result but consider that the initial state is it- 
self the SP layer, with a thickness that is proportional to the 



square of the resistivity, 



It immediately follows 



that the fastest growing mode of the secondary island insta- 
bility scales as k oc ?7~ 3 / 8 and the corresponding growth rate 
as 7 oc ?7 _1 / 4 . These simple estimates agree remarkably with 
the more in depth analysis by Loureiro et al. (2007|. 

The argument can be carried further in the non-linear 
regime to estimate the non-linear reconnection rate during 
the subsequent growth of the islands. In the case of the tear- 
ing mode, the island growth is governed by the Rutherford 
equation (Rutherford 1973| l, stating that the size of the mag- 
netic islands (wj) grows in time linearly (rather than expo- 
nentially). The temporal derivative of the size of a magnetic 



island is wi oc rjA' / no ( |Kadomtsev 1992 p. 76), where, for 
the initial configuration assumed above, 



A'= 2 (^-ka 
a \ka 



(8) 
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At long wavelengths, where ka << 1, A' oc a~ 2 , and assum- 
ing again the initial thickness to be that of the SP layer, the 
Rutherford equation leads to a rate of island growth indepen- 
dent of resistivity, u)j oc rp. 

Note that the dependence on k is not considered in the 
scaling argument because the secondary islands are observed 
in the simulation to coalesce and the value of k is expected to 
settle on a value constrained more by the macroscopic system 
rather than by the fastest growing linear mode. 



5 Role of the flow pattern in unsteady reconnection 

A consequence of the onset of the secondary island instabil- 
ity is the profound modification of the plasma flow patterns 
in the reconnection region. To investigate the flow, a specific 
simulation is used. The initial state is a Harris sheet with 
magnetic field given by eq. (|5j, uniform initial temperature 
corresponding to v s /va = 1 and density 



but the same series of events develops also with open bound- 



p = p /cosh 2 (z/a)+p b 



(9) 



with pi, — po/10. The evolution is followed with visco- 



resistive MHD, using the FLIP3D-MHD code ( jBrackbiU 
[T99T) . 

The simulations have uniform resistivity corresponding to 
a Lundquist number, S = tjoCLVa/t] = 10 4 and a viscosity 
corresponding to a Reynolds number, R — poVAd/rj = 10 4 . 
All results are presented in normalised units where the mag- 
netic field is normalised to Bq, the density to po, time with 
the Alfven time ta — ci/va, space with the initial current 
sheet thickness a and velocity with the Alfven speed va = 
Bof(poPo) 1 ^ 2 - The simulation box has sides L x /a = 240 
and L z = 60. The code is 3D but the y dependence is sup- 
pressed. 

During the SP regime, the flow is primarily directed nor- 
mal to the sides of the SP layer. In the two inflow sides, 
the flow is vertical, and in the two outflow sides it becomes 
horizontal, veering within the SP layer. Figure [5]-a shows 
such a configuration for the present simulation. The actual 
flow in the simulation is just as the typical textbook SP car- 



toon (Biskamp 2000). 

As the simulation evolves, the flow pattern changes and 
becomes primarily focused near the separatrices. This is a 
typical feature of the flow in presence of magnetic islands 
due to the tearing instability. Figure |5]-b and Fig.|5]-c demon- 
strate the transition. A x-point is forming in the center as 
a consequence of the secondary island instability. The flow 
pattern becomes now just like the textbook cartoon of the 
flow near an x-point in the tearing mode (Biskamp, 2000). 

Still later times are shown in Fig. [6] After the onset of the 
secondary island instability progressively more secondary is- 
lands and x-points are formed. The new formed islands coa- 
lesce with pre-existing islands and migrate with the flow to- 



ary conditions (Skender and Lapenta 2010 Bettarini and 
Lapenta 2010 ) with the only difference that the flow carries 



the formed islands outside the open boundaries. 

The presence of multiple islands causes the flow to form 
characteristics loops. In a previous paper (Lapenta 2008), 



the name conveyor belts was used to explain the role of such 
flow loops. The flow assumes a pattern where the flow lines 
that come out of the reconnection process bend and return 
towards the same reconnection region. The effect is sim- 
ilar to that of recursive reconnection (Par nell et al.| 2008) 
observed in magnetic skeleton configurations representative 
of solar coronal configurations. In recursive reconnection 
the outflow of one reconnection site feeds into another. In 
the present simulation the same recursive process is also ob- 
served. In Fig [6[c the outflow from the x-point located at 
x/a = 115 feeds the reconnection process at the x-point lo- 
cated at x/a = 90 (and the similarly by symmetry on the 
other side). But more conspicuously, the same reconnection 
x-point at x/a = 115 feeds itself with most flow lines going 
out of it and coming back forming a conveyor belt. 

In accordance with the theoretical predictions above for 
reconnection to be fast the flow must become stronger in the 
unsteady regime. Comparing the scales in Fig. [5] with Fig. [6] 
the flow accelerates strongly. Figure [7] shows separately the 
inflow component along z and the outflow component along 
x. The outflow is at the Alfven speed and the inflow is a 
significant fraction of it. At steady state the inflow speed 
equals the reconnection electric field and is by definition the 
reconnection rate. 

The concept of streamlines and in particular of conveyor 
belts should not be confused with actual recirculation of 
plasma. The flow pattern is distinct from the matter actu- 
ally flowing. On those stepping into rivers staying the same 



wards the two ends of the simulation box (Lapenta 2008 1. 
The present simulation uses periodic boundary conditions 



other and other waters flow ( |Heraclitu s of EphesusJ |2001 
DK22B12,Fragment 91). The flow can present the same 
pattern over a certain period but the plasma elements going 
through the flow are other and other. The streamlines are in- 
stant photos of the velocity field at a given instant. They are 
obtained here with the paraview software using the fourth- 
fifth order Runge-Kutta ODE solver that traces the lines ev- 
erywhere tangent to the velocity field at that instant. But an 
actual fluid element feels at every instant the local speed at 
that instant. A moment later it feels another speed. Flow 
lines are not the same thing as particle trajectories. To anal- 
yse actual trajectories of fluid elements we use Lagrangian 
markers. Lagrangian markers are real particles of the fluid 
that move with the local instantaneous velocity just as a phys- 
ical fluid would. These are not kinetic particles moving with 
their peculiar speed but rather fluid elements moving with 
the local flow speed. The FLIP3D-MHD code is based on 
a dual grid-particle representation that includes the presence 
of Lagrangian markers ( |BrackbilT] 1 1 99 1 ) ) . 

Figure [8] shows a selected number of trajectories of La- 
grangian markers. The initial starting point at time t = 
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in the simulation is chosen appropriately to identify parti- 
cles that will go through multiple reconnection sites. Some 
trajectories only pass reconnection once: those are particles 
that are drawn toward the neutral line (z /a — 30) and then 
ejected away from it (reddish trajectories on the left and light 
cyan on the right). However, there are other classes of tra- 
jectories (green on the left and magenta on the right) that are 
first drawn towards the neutral line and after being pushed 
away, come back a second time towards the neutral line and 
then are finally ejected away from it. These are particles that 
are captured in the conveyor belts and are pushed towards 
a reconnection site and then away from it and then towards 
another one again. Note, however, that these trajectories are 
still monotonic in the horizontal direction: the same parti- 
cle does not turn back to go towards the same x-point. For 
example on the left of center, the markers first take part to 
reconnection at the site located near x/a = 115 and then at 
x/a = 90. No particle is detected to go through the same 
reconnection site twice. Just as the citation by Heraclitus 
suggested. 

6 Discussion 

A magnetic field embedded in a laminar perfectly conduct- 
ing fluid preserves its topology for all time (Parker 79). Al- 
though ionized astrophysical objects, like stars and galactic 
disks, are almost perfectly conducting, they show indications 
of changes in topology, "magnetic reconnection", on dynam- 
ical time scales (Parker 1970, Lovelace 1976, Priest & Forbes 
2002). Reconnection can be observed directly in the solar 
corona ( Innes et al 1997, Yokoyama & Shibata 1995, Ma- 
suda et al. 1994), but can also be inferred from the existence 
of large scale dynamo activity inside stellar interiors (Parker 
1993, Ossendrijver 2003). Solar flares (Sturrock 1966) and 
7-ray busts (Fox et al. 2005, Galama et al. 1998) are usually 
associated with magnetic reconnection. Previous work has 
concentrated on showing how reconnection can be rapid in 
plasmas with very small collisional rates (Shay et al. 1998, 
Drake 2001, Drake et al. 2006, Daughton et al. 2006), which 
substantially constrains astrophysical applications of the cor- 
responding reconnection model^] 

In comparison, our results point out to an attractive uni- 
versal astrophysical solution for the magnetic reconnection: 
magnetic reconnection in presence of turbulence is always 
fast independently of the collisional or collisionless state of 
plasmas. The formal criterion of the applicability of theLV99 

2 We note that if magnetic reconnection is slow in some astro- 
physical environments, this automatically means that the results of 
present day numerical simulations in which the reconnection is in- 
evitably fast due to numerical diffusivity do not correctly represent 
magnetic field dynamics in these environments. If, for instance, the 
reconnection were slow in collisional media this would entail the 
conclusion that the entire crop of interstellar, protostellar and stellar 
MHD calculations would be astrophysically irrelevant. 



model based on the MHD approximation to realistic plasmas 
is discussed in Eyink, Lazarian & Vishniac (201 1). 

Our paper provides a connection between the LV99 the- 
ory and numerical simulations demonstrating the sponta- 
neous onset of fast reconnection. Indeed, according the LV99 
model the reconnection in the initially laminar state is slow. 
The outflow, however, incites turbulence in the system. This 
turbulence induces faster reconnection and a wider outflow 
gets more turbulent. This should induce a positive feedback 
process which results in the flare of reconnection. We be- 
lieve that this is what is the actual theoretical explanation of 
the spontaneous onset of fast reconnection observed in pa- 
pers by Lapenta et al. 

The spontaneous onset of reconnection is likely to be the 
cause of the solar flares. The flares require the accumulation 
of magnetic flux of opposite polarity prior to the flare. The 
accumulation means low reconnection rate, which is readily 
achived for the observed mean level of turbulence in solar at- 
mosphere (Eyink et al. 201 1). As the flux is accumulated the 
effects that we described in this paper are expected to come 
into place and a abrupt release of energy through magnetic 
reconnection is expected. A similar process can take place 
during gamma ray busts (Lazarian et al. 2002, Zhang & Yan 
2011). A relatively low level of turbulence is a prerequisit 
for the flares. In the presence of high amplitude turbulence, 
reconnection is going to be fast from the very beginning pre- 
venting the accumulation of magnetic flux at the pre-flare 
stage. 

Observationally, predictions of the LV99 model for flares 
are in rough agreement with the estimates of the thickness 
of the outflow regions obtained in observations (Ciravella & 
Raymond 2008). In addition, the initiation of the reconnec- 
tion by Alfvenic perturbations predicted in LV99 was con- 
firmed by observations in Sych et al. (2009). Direct observa- 
tions of turbulent reconnection in other environments, e.g. in 
interstellar medium are not feasible at the moment. However, 
the indirect evidence of LV99-type reconnection comes from 
the observations of magnetic fields in molecular clouds and 
accretion disks. In fact, on the basis of LV99 model Lazar- 
ian (2005) predicted the process of magnetic field removal 
termed "reconnection diffusion". Numerical studies of the 
reconnection diffusion (see Santo de Lima et al. 2010, 201 1) 
deliver results on expected magnetic flux diffusion that are 
consistent with observations. 

The caveat here is that the LV99 model is the 3D model 
and it appeals to the 3D wondering of magnetic field lines. 
We believe that the fact that the acceleration was observed 
also in 2D systems testify that the modification of LV99 
model is applicable to the systems of reduced dimentionality. 
It is discussed in Eyink et al. (201 1) that the existence of the 
"rough" turbulent velocity field should induce the Richarson 
diffusion and therefore the effective magnetic field wander- 
ing. The corresponding theory have not been developed for 
the 2D turbulence. However, results on the acceleration of re- 
connection in the presence of turbulence (Loreiro Kulpa- 
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Dubel et al. 2010) suggest that the results of LV99 may carry 
over to 2D systems. However, as turbulence in 2D and 3D is 
different, we feel that it does not worth rederiving the LV99 
results for 2D, but it is better to concentrate on matching the 
quantiative results for the actual real world 3D reconnection 
as numerical simulations are getting sufficiently powerful to 
study the effects in 3D. This is the theme of our intended next 
paper. 

Our paper sends a warning signal to a naive interpretation 
of astrophysical reconnection within laminar reconnection 
models. We showed that pre-existing turbulence as well as 
turbulence induced by the process of magnetic reconnection 
itself can dramatically change the reconnection rates. Inter- 
estingly enough, the small - scale reconnection events may 
be determined by small scale physics, but the total reconnec- 
tion rates may still be determined by turbulent level in the 
system. 

7 Summary 

The results of the present paper can be briefly summarised 
in the following way: 

1. The LV99 theory of fast turbulent reconnection and 
more recent results on the spontaneous onset of turbulent re- 
connection are related to each other. We claim that the flares 
of reconnection observed in simulations are due to the feed- 
back by turbulence induced by the outflow on the reconnec- 
tion process. 

2. The numerical results suggest that flares of magnetic 
reconnection, e.g. explaining solar flares, may be explained 
on the basis of pure resistive MHD. 

3. Our work testifies in favor of a widely applicable model 
of astrophysical reconnection that is independent of subtle 
plasma effects. At the same time it send a warning to at- 
tempts to model astrophysical reconnection without includ- 
ing effects of turbulence. 
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Fig. 5. Early evolution. The flow lines are shown superimposed over a false colour representation of the flow speed. The SP layer first forms 
and the flow pattern around it changes its nature. Three times are shown. Blow up of the central region, the full box is 240x60. 



Lapenta and Lazarian: Fast Turbulent Reconnection 



13 



35 

30+ 

25 



a) v a t = 332 



U Magnitude 
0.04 0.08 0.1 0.20.2 



0.0003 



0.22 




b) v a t = 352 



U Magnitude 
0.04 0.08 0.1 0.20.2 



0.0003 



0.22 




c) v a t = 368 



U Magnitude 
0.2 0.4 



le-6 



0.53 




13 - VI 



: 

■ — — MoU ■ 1 QU 



Fig. 6. Later evolution. The flow lines are shown superimposed over a false colour representation of the flow speed. The SP layer is 
destabilised by the secondary island instability. Three times are shown. The first two have the same color scale, but the third has higher 
values to avoid saturation. Blow up of the central region, the full box is 240x60. 
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Fig. 7. Flow pattern at time t/TA = 375.5. The outflow component U x is shown in panel a and the inflow component z in panel b. Blow up 
of the central region, the full box is 240x60. 
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Fig. 8. Trajectory of Lagrangian markers. Miltiple trajectories are 
shown in different colours. Each trajectory starts at a circle and ends 
at a cross. The initial positions for the aprticles are equally s paced 
between x — L x — ±7L and x — L x /2 = ±14L on either side from 
the center of the x axis and at vertical position z/L = 29.9167, just 
a little below the central neutral line. 



